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ABSTRACT 


This report describes results of a NASA sponsored study to investigate 
and analyze the general aviation air traffic flow patterns at uncontrolled 
airports and to develop traffic pattern concepts which minimize the mid-air 
collision hazard in uncontrolled airspace. An analytical approach to 
evaluate mid-air collision hazard probability as a function of traffic 
densities is developed which is basically independent of path structure. 

Two methods of generating space-time interrelationships between terminal 
area aircraft are presented. One is a deterministic model to generate pseudo- 
random aircraft tracks and this is compared with analyses of results from 
available real data. Some hazard measures are presented for selected traffic 
densities, A second model which is a statistical model in preliminary form 
is discussed. An analytical expression for the statistical description of 
procedure path deviations is derived. 

An analysis of a developed generalized expression for mid-air-collision 
(MAC) probability results in the conclusion that the probability of encounter- 
ing a hazard should be minimized independently of any other considerations 
and that given a certain encounter rate the number of encounters involving 
visible-avoidable aircraft should be maximized at the expense of encounters 
in other categories. It is shown that the pilot "look-time" should be pro- 
portioned according to the magnitude of the probability of finding a visible- 
avoidable hazard. 
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CHAPTER 1 


INTRODUCTION 

In the growing air traffic environment of today the investigation of 
mid-air collision hazard is of prime significance to our society. Of 
particular interest is the relationship between prescribed procedures and 
the mid-air collision hazard to general aviation aircraft In uncontrolled 
terminal airspace. With the objective of developing traffic pattern concepts 
to minimize mid-air collision hazard In uncontrolled airspace, NASA-Wallops 
personnel have developed a data base consisting of position-time tracking 
data. The data provide information on terminal area air traffic patterns 
including those flown by general aviation aircraft at several medium-sized 
city airports thal generally represent uncontrolled airport environments. 

The Research Triangle Institute, under contract NAS6-2312, is providing 
support by developing an analytical concept for the analysis of data obtained 
by NASA. Various approaches have been invest i r at<»d to provide realistic 
measures of mid-air collision risk and reliable methods for mathematical 
comparison of different air traffic flow situations. These have included the 
development of preliminary computer mathematical models of the uncontrolled 
terminal area traffic and methods for validation of these analytical models 
using the empirical air traffic data collected by NASA-Wallops personnel. 

1.1 Problem Statement and Significant Issues 

Published procedures for uncontrolled terminal area environments pre- 
scribe that an aircraft upon approaching the runway must generally enter a 
rectangular racetrack-type pattern which involves left turns at relatively 
high bank angles in a counter-clockwise runway encircling path. Pattern 
altitudes may vary from runway to runway but generally gradually descending 
legs are prescribed to allow orderly progression to touchdown. This basic 
procedure, illustrated in Figure 1-1 (see ref. 1), provides tlie pilot 
under VFR conditions with a good continuous view of the runway during the 
approach and landing phase. In an uncontrolled environment each pilot is 
dependent basically on the "see-and-avold' 1 concept to prevent mid-air 
collisions so that the question of whether or not prescribed terminal 
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Approach Pattern 




approach procedures are consistent with minimizing mid-air collision hazard 
in a see-and-avoid environment is fundamental. 

To answer thlB question It Is necessary to evaluate the present un- 
controlled patterns flown and to determine Improvements in a pilot's 
ability to see another aircraft (assuming he looks) for various changes in 
the traffic pattern concept. In refs. 2 and 3 some analysis of approach 
data has Indicated typical adherence to procedure in an uncontrolled 
environment. Additionally consideration is given to pilot visibility from 
the cockpit through the use of aircraft vision envelopes to give insight 
into how approach pattern visibility might be evaluated. No evidence has 
been found of previous studies which evuluate the approach pattern on the 
basis of mid-air collision risk in a see-and-avoid environment. However, 
this type procedure has been proposed as a technique for evaluating the 
effect of VFR towers on. traffic flow and sequencing (ref. 4), 

In a controllvd ur.vii orment separation standards (or minimums) are 
normally enforced n that '’natural" course deviations will allow very 
little collision ri.-k (see refs. 5 and 6). The idea of collision frequency 
has been developed (ref. 5) as the frequency with which the difference 
between the actual positions of two aircraft comes within a certain value. 

A prediction of this frequency can then be termed a probability of collision 
The task of relating collision risk to a traffic configuration can be broken 
into two parts: (1) determining the frequency with which aircraft are 

exposed to risk by passing close together and (2) determining what chance 
of collision is inherent i,i the passing. One approach to evaluate a given 
path structure is to analyze the flight paths and calculate the relative 
frequency that these paths are "close" considering total flight time and 
assuming that aircraft proceed as if no collision occurred when path cros- 
sings do occur. In ref. 5 this is called "blind-flying collision r.sk." 

In the uncontrolled region, collision hazard to a particular aircraft 
has been stated to depend on (1) density of aircraft, (2) heading distribution 
and (3) speed distribution (ref. 5). For the situation involving random 
independent paths a small deviation in flight path will not change the 
environment and therefore not change the hazard. The statement has been 
made that the chance of having a collision is then determined by the 
environment and is unaffected by deviations from the flight path. For the 
situation of interest here, the assumption of random independent paths is 
not valid since there is a prescribed procedure even in the uncontrolled 
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terminal area. Preliminary Indications (ref. 2) are that prescribed path 
deviations are not small enough to be negligible and in fact can contribute 
to the mid-air collision risk. 

In ref. 7 collision hazard is actually divided into two parts. A 
measure of pilot workload is consiJ .red to be related to the percent of 
time that a given aircraft finds another aircraft in an "encounter" status. 

An encounter status can be defined in a number of ways and can include re- 
lative range, closing velocity, and closing acceleration. A second contri- 
butor to collision hazard is encounter rate, which can be related to pro- 
bability of an encounter during some given time period in which this time 
period is very small compared to the 4,ime basis for measurement of the 
encounter rate. 

In evalu~i-1 \g the inherent risk ro be associated with a particular 
traffic pattern, it seems apparent that the important factors for considera- 
tion should include (1) density of aircraft in the terminal area, (2) conflict 
points inherent in the traffic patterns, (3) speed distribution and type of 
aircraft, (4) visibility, and (5) pilot workload. These factors are con- 
sidered in developing an analysis of the uncontrolled terminal area environ- 
ment in terms of a collision hazard probability. 

1.2 Summary of Results 

A generalized expression for mid-air-collision (MAC) probability is 
developed which yields two basic hypotheses about the structuring oi patterns 
for the terminal area. First, the probability of encountering a hazard 
should be minimized independently of any other considerations. Second, and 
perhaps not as obvious, if given a certain encounter rate, the number of 
encounters involving visible-avoidable aircraft should be maximized, thus 
minimizing the number oi encounters in other categories (invisible-avoidable, 
invisible-unavoidable, or visible-unavoidable). Additionally, an Important 
consideration concerning pilot function is borne out. The pilot "look time" 
should be proportioned according to the magnitude of the probability of 
finding a visible-avoidable hazard. Therefore traffic pattern design should 
(1) minimize encounter rates and (2) maximize visible-avoidable encounters 
under the constraints of aircraft cockpit binocular design. 
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Two approaches have been Investigated and are described in this report 
relating to the required analysis of air traffic patterns and their relation- 
ship to hazard in ti.e uncontrolled terminal area environment. An analytical 
occupancy model based on a statistical description of the terminal area 
airspace is presented. This can serve as a basis for evaluating various 
approach patterns once the procedure is completely developed and evaluated. 
Particular pattern characteristics can be defined in terms of input statistics 
and a measure of hazard probability can be analytically determined. The 
resulting statistical measure of hazard probability can be evaluated for 
various pattern characteristics. 

Secondly, a trajectory model is described which involves a set of para- 
meters which can be varied to generate approach paths similar to those 
present in real data. Thin deterministic approach is demonstrative of a 
capability to generate realistic traffic paths using a computer algorithm 
which can readily be altered to conform to essentially any desired approach 
pattern. 

The deterministic model is used in a preliminary term to demonstrate 
evaluation of encounter probabilities for a path structure designed o 
represent that currently used in the uncontrolled environment. Paths are 
time correlated by assignment of touchdown times to entering aircraft 
according to a Poisson distribution with a rate parameter dependent on 
the assumed aircraft density. This can provide a baseline measure of the 
hazard associated with a particular traffic pattern. Flight path deviations 
are measured and distributions of across-path errors are compared with 
those obtained from the Wallops Island data analysis to support the premise 
that this modelling approach can be used to generate realistic flight paths 
representing the uncontrolled terminal area environment. 

Finally, recommendations to improve the deterministic model are given 
along with suggested extensions to the analytical occupancy model approach. 
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CHAPTER 2 


MID-AIR COLLISION MODELS 

The ultimate value of a hazard model for uncontrolled airspaces lies 
in the utility of the model for predicting and minimizing the mid-air 
collision probability. In this chapter a general expression for the 
mid-air collision probability is advanced. The expression shows explicit 
dependence on aircraft visibility and detectability, pilot function, and 
hazard probability. The expression Illustrates rather straightforward 
minimization tactics that can be pursued, and dramatizes the dominance of 
the hazard probability. Models for hazard probability are developed in 
Chapter 4. 


2.1 Mid-Air Collision Probability 


An abstract aircraft space is depicted in Figure 2-1. The space 
dichotomizes aircraft into those that are non-hazardous (H) to a target 
aircraft and those that are hazardous (H) and would result in a mid-air 
collision if no evasive action were taken. The subset of hazardous air- 
craft has been further refined into disjoint subsets H_^ corresponding to 
generalized geometrical locations (including space and time, for example) 
of hazardous aircraft. It is convenient to denote tha subset of hazardous 
aircraft by 


H = 


n 

£ H 

i=l 1 


( 2 - 1 ) 


where n is the number of disjoint subsets H^. Now let MAC denote the 
subset of hazardous aircraft that result in mid-air collisions. Since 
the set MAC is a subset of H, and H is a finite union of disjoint sets, 
this set caii be written as 

n 

MAC - MAC OH, . (2-2) 

i=l 1 




7 



Abstract Aircraft 
Space 


mm 




_i Subspace of Nonhazardous Aircraft 
HH Subspace of Hazardous Aircraft 

lj | Subspace of Hazardous Aircraft 
Resulting in Mid-Air Collisions 


Figure 2-1. Aircraft Space. 


8 




The probability of a raid-air collision Is 


n 

PfMAC] = ^ 

i-1 

P [MAC OH.] 

n 

- E 

i«l 

P[MAC/H ± ] P[H ± ] 


where P[H^] is a hazard probability. Denoting the complement of MAC by 
MAC, note that P[MAC] can be written 


P [MAC] = 


n 


E 

-1=1 


(1 - PfMAC/H.jD Pfiy 


n 

■ P[H] - P[MAC/H i ] P£H ± ] . (2-4) 

1=1 


In words, eq. (2-4) says the probability of a mid-alr collision is the 
probability of encountering a hazard, P[H], minus the probability of 
encountering and avoiding the hazard. If P[MAC/H^] = 1 for all i (i.e., 
all hazard conditions are avoidable with probability 1) then P[MAC] = 0. 

A certain subset of the set H may contain aircraft that are 

invisible, unavoidable, or both. Such hazardous aircraft result in mid-alr 

collisions with probability one. Thus eq. (2-4) can be rewritten 


P[MAC] = P[H) - ^ PIMAC/^J P[H ± ] (2-5) 

iel 

where H' = H. denotes the set of hazardous aircraft that are visible (V) 
iel i 

and avoidable (A) . 


9 



The set (or event) MAC is the event that a visible, avoidable, 
hazardous aircraft is "looked at" (L) by the pilot (but not necessarily 
seen if the hazard is not perfectly detectable) and detected (D) . It is 
assumed that a visible, avoidable target that is detected will be avoided, 
thereby ruling out of this analysis the Irresponsible pilot who fails to 
take evasive action. The notation for MAC is 

MAC - LHd . (2-6) 

Equation (2-5) can now be written 

P [MAC] = P[H] - 2 PtLODOH^ 

isl 


- P [H] - 2 PID/LOHj,] P[L/H ± ] PD^] 

iel 

- P[H] - ^ P [L/H ± ] PtHj,] • (2-7) 

iel 

The last simplification follows from the observation that the inherent 
detectability of a hazardous aircraft is independent of whether or not the 
pilot looks at it. The interpretation of eq, (2-7) is the same as the 
interpretation following eq. (2-4). However the additional interpretation 
can be made that the probability of encountering and avoiding a hazard 
is given by the probability of encountering, looking at, and detecting a 
visible, avoidable hazard. Several plausible minimization tactics follow. 

2,2 Minimization Tactics 

It seems plausible that the summation term in eq. (2-7) should be 

independent of (or weakly dependent upon) the probability of encountering 

a hazard, P[H]. That is, the nature of the set H' = H. of visible and 

iel 1 

avoidable hazards, and the probability of encountering, looking at, and 
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detecting this set of hazards, should be independent of P[H]. With this 
assumption it can be concluded that P[H] , the probability of encountering 
a hazard, should be minimized independently of any other considerations. 

This is of course a traffic pattern and scheduling issue. Next, the set 
H* of visible, avoidable, hazardous aircraft should be maximized, subject 
to the consideration that hazards which are likely to be encountered, looked 
at, and detected (corresponding to large values of P[D/H^] PfL/H^] P[H^]) are 
to be included at the expense of others. If P[D/H^] P[L/H^] P[H^] is inde- 
pendent of H^, then the cardinality of the set of visible, avoidable targets 
should be maximized independently of other considerations. 

Finally, the following interesting minimization with respect to pilot 
function is evident. Let g(i) = P[D/H^] P [ ] denote the probability of 
encountering and detecting the visible, avoidable hazard H^. Then to 
minimize P[MAC], maximize the objective function 


J = £ g(l) PEL/Hj,] (2-8) 

iel 


where P[L/H^] is the probability of looking at hazard This function is 

maximized by choosing (Schwartz inequality) 


P[L/H ± ] = k g(i) (2-9) 

where k ^ g(i) (assuming the pilot looks in some generalized 

hazard location with probability 1). Thus the pilot should spend [g (1) / 
g(i)]% of his time looking in hazard location H^. Roughly speaking, the 
percentage of time spent looking in hazard location should be proportional 
to the probability of detecting a visible, avoidable hazard in that location. 
The probability of detection of a visible, avoidable hazard is maximized 
through landing pattern design and aircraft binocular design (ref. 3), 
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CHAPTER 3 


BASIC AIRSPACE MODEL 

The critical issue in a model for mid-air collision probability 
seems to be the determination of hazard probability, PfE^]. Several hazard 
probability models have been proposed in relation to collision avoidance 
studies (see, ref. 8 and the references therein). These models are 
generally geometrical in nature, as they must be, involving complicated 
kinematic relationships between two or more aircraft. To date, the 
models have not Included the effects of random aircraft entrance times into 
controlled or uncontrolled airspaces, nor have they accounted for random 
aircraft trajectories within an airspace. (The suggestions contained in 
ref. 9 represent a notable exception to this statement.) 

In this chapter the issue of random aircraft entrances into uncontrolled 
airspaces is addressed by developing a basic airspace model. First, consider 
a simple, but representative, landing pattern and then pose a homogeneous 
Poisson model for the number of aircraft that land in an arbitrary time 
Interval [o,t). By assigning aircraft entrance probabilities for upwind, 
downwind, and crosswind legs of the pattern (based on published data, e.g„, 
ref. 8), a basic airspace model for the number of aircraft occupying any 
region of the pattern at any point in time can be obtained. The airspace 
model is then generalized to Include the nonhomogeneous effects of daily 
fluctuations in traffic density. Some comments are advanced concerning 
extensions of the basic airspace model to Include random three-dimensional 
aircraft trajectories. 

In Chapter 4 several hazard probability models are derived from this 
basic airspace model. In Chapters 5 and 6 the airspace model is used to 
partially specify a trajectory simulation program for generating random 
aircraft trajectories and obtaining various statistical measures of an 
uncontrolled airspace. 


3.1 Homogeneous Poisson Model for Aircraft Landings 


Let the number of aircraft, N(t), landing in the time interval [o,t) 
be a homogeneous Poisson counting process with intensity v and rate parameter 
X = vtj i.e. , 


PRECEDING PAGE BEBNT 
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Aircraft Entering 
/ at Leg 3 

Travel Time: T 2 

Velocity: V 2 (t) 


Figure 3-1. Representative Landing Pattern. 
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PlN(t) - n] 


(3-1) 


(vt) n -vt 
e 


where P[N(t) ■ n] denotes the probability that N(t) « n and vc Is the 

average number of aircraft landing in [o,t). 

Now consider Figure 3-1 depicting a five-leg landing pattern. One 

might reasonably ask what the probability is that N m (t) aircraft will enter 

the landing pattern at leg m in the time interval [o,t). To answer this 

question, assume the general landing pattern contains M legs, aircraft 

enter leg m with probability p^, and p^ - 1. Then the probability 

that N (t) « n can be written ,i! 

m m 


P[N (t) - n ] 
m m 


E P[N(t) * n/N(t + t ) - N(t ) - n] 
m m m m 


• P[N(t + t ) - N(t m ) = n] 


(3-2) 


where T. is the appropriate "travel time" required for an aircraft to 
traverse the j leg, and t m = J'' T^ is the total time required for 
aircraft to land after entering-^ t{?e landing pattern at leg m. The 


an 


assumption here is that every aircraft that enters leg m in the interval 

[o,t) must land in the Interval [t , t + t ), where t is the time required 

for an aircraft to travel from the entrance of leg m to the landing point. 

Then eq. (3-2) follows from the law of total probability. 

The first term in eq. (3-2) is the probability of n^ successes (a 

success defined here as an entrance at leg m) in n trials (a trial defined 

here as a landing in the time interval [t , t + t ). The distribution is 

m m 

b inomial : 
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(3-3) 


P[N(t) - n/N(t + t ) - N(t )] « \ \ (1 - pj n “ n m . 

m m tn m \n/p 

\ m/ r m 

The second probability statement is Poisson with rate parameter X(t + t m ) - 
Xt m , by virtue of the fact that increments of a Poisson process :,*e also 
Poisson: 

P[N(t + t ) - N(t )] - - X ^ ~ e" Xt . (3-4) 

m ra n I v ' 

Equation (3-2) is now evaluated as follows: 


P[N (t) 
m 


n ] 
m 


V _2JL 

(n-n ) In t 
n=n mm 

m 


n 

ID * ~ 

(1 ' 


x n-n Xt n -Xt 
p ) m — r - ~ - e 


nl 


n 


m 


a - P j J 


m /,.vn -Xt \ 1 r m' ,,^1 

*= - — (Xt) m e / , 


n 


ml 


1-0 


11 


(Xt) J 


(p Xt) n m , 

*m -P Xt 

e 


n 


(3-5) 


m 


The conclusion is that the number of aircraft entering leg m is Poisson 
with rate parameter p^ Xt: 

N (t) : P[p m Xt] . (3-6) 

tn m 

The parameter p m can be obtained from experimental observations regarding 
the relative frequency of aircraft entrances at leg m (ref, 3). The para- 
meter X is also obtained from experimental observations. 
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These calculations are based on preservation of the property that the 

average number of aircraft landing in the interval [t^, t + tj) is the 

average number entering leg 1 in the interval [o,t) plus the average 

number entering leg 2 in the interval [t^ - t^, t + tj- t 2 ) ••• plus the 

average number entering leg M in the interval [t^ - 

Since the average number of aircraft entering leg m in the interval 

[t, - t , t + t. - t ) is p Xt, it follows that 

1 1 nr 1 m r m ’ 


V + t l ~ V* 


M 


m*l 


p Xt 
m 


Xt 


(3-7) 


and the averages are preserved. 

Actually, more can be shown. Consider the Individual Poisson 
processes N^. The number of aircraft landing in [t^, t + t^) 

equals the number of aircraft entering leg 1 in the interval [o,t) plus the 
number entering leg 2 in the interval [t^ - t 2> 11 + ^ plus the 

number entering leg M in the interval [t^ - t^, t + t^ - t^) ; 

N(t + t L ) - NC^) * N^t) + [N 2 (t + t L - t 2 ) - N(t x - t 2 )J 

+ ... + [N M (t + - t M ) - N(t 1 - t M )] • (3-8) 

The characteristic function of N(t + t^) - N(t^) is 


*(«) = E[e « t+t l> 


e ^ e -wt (Xt) u e jwn 
e Xt(e^ w - 1) 


(3-9) 
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But this can ba written 


M 

2 P m " 1) * 

. / \ m*l 

Kw) ■ e 

M o Xtfe*^ - i/ 

- H e Pm A ^ ' (3-10) 

m*l 

The ui term of the product is the characteristic function of the Poisson 

random variable N (t + t- - t ) - N (fc, - t m ) . Thus it follows that the 
m 1 m m 1 m 

N^(t), N^(t) are independent Poisson random processes. These results 

will be generalized for nonhomogeneous airspaces shortly. 

The conclusion is that Poisson models for the independently arriving 
aircraft at the various landing legs lead to a Poisson model for the 
number of aircraft landings on an arbitrary time interval (and vice-versa), 
provided relevant time intervals are carefully chosen. The mathematical 
model is therefore consistent with the physical constraint that aircraft 
which enter the pattern must land after an appropriate delay. The appro- 
priate delay has been assumed deterministic, corresponding to the assump- 
tion that aircraft trajectories are deterministic. It Is very difficult 
to relax this requirement and still obtain tractable models, but in 
Chapter 4 some preliminary results are advanced for random aircraft 
trajectories. 


3.2 Interval Occupancy in Homogenous Air Spaces 

An important determination to be made for the calculation of hazard 
probabilities is the number of aircraft, Q(x;t), located in the spatial 
Interval [x, x + Ax) along the landing pattern at time t (see Figure 3-1). 
Consider an interval Tx, x + Ax) on the first leg (i.e., d^ £ x < x + Ax £ d£, 
with d^ the beginning of leg 1 and the beginning of leg 2) . It follows 

that Q(xjt) equals N^(t - t^(x + Ax)) - N^(t - t^(x)), the ' ..mber of aircraft 
entering the landing pattern at leg 1 in an appropriate interval [t - 
t^ (x + Ax), t - t^Cx)), The time C^x) is the time required for an aircraft 
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with velocity v^(t) to traverse the distance from to x: 



v, (t) dt - x - d 


(3-11) 


Thus Q(x;t), for d^ £ x x + Ax _< d^, Is Poisson with rate parameter 
P 1 vft-^x + Ax) - t^x)]: 


Q(xj t) : P[p^v(t^(x + A* 1 :) > t^(x))] , d^ £ x <_ x + Ax <_ • (3-12) 


When the velocity is a constant v , then t, (x) ■ (x - d_)/v and the rate 

ox x o 

parameter becomes p 1 vAx/v q . The probability that Q(x;t) » q is 


(p.vAx/v ) q -p.vAx/v 

P[Q(x;t) “ q] = j e , d^ <_ x < x + Ax <_ (3-13) 


For an Interval [x, x + Ax) on the second leg, Q(x;t) is the number of 
aircraft entering the landing pattern at leg 2 in an appropriate interval 
[t - t^(x + Ax), t - t 2 (x)), plus the number of aircraft that entered leg 1 
in an appropriate interval [t - t^(x + Ax), t - t^(x)) and continued in the 
pattern from leg 1 to leg 2. The time t 2 (x) is the time required for an 
aircraft with velocity v 2 (t) to traverse the distance from d 2 to x > d 2 ; 

t 2 (x) 

J ' v 2 (t) dr = x - d 2 . (3-14) 

o 
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The time t 1 (x) is simply t-^x) - t 2 (x) + T.^ where is the time required 
for an aircraft to traverse leg 1. Thus 

Q(x;t) - [N 2 (t - t 2 (x)) - N 2 (t - t 2 (x + Ax))] 

+ [N 1 (t - t 2 (x) - Tj) - N x ( t - t x (x + Ax) - T x )] . (3-15) 


By the independence and homogeneity of and N 2 , it follows that Q(x;t) is 
Poisson with rate parameter (p.^ + p 2 ) \) [t 2 (x + Ax) “ i,e, » 


Q(x;t) ; PKp^ + P 2 ) v ( t 2 ^ x + Ax) “ t 2 (x))] (3-16) 


When v^(t) =• v 2 (t) = v q , 


Q(xjt) 



(3-17) 


Continuing in this manner, it is easy to show that for x rontained in 
th 

the m leg, d <x<x+Ax<d , 
m — — ra + 1 


Q(x ; t) 



p,v(t (x + Ax) 



(3-18) 


For v 1 (r) = v 2 (x) v ra (x) - v q , 
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Q(x;t) : P 


(3-19) 


£ Pj V 4*/v„ . 

J-l 


m 

The average number of aircraft in [x, x + Ax) is v(Ax/v q ) 23 p ^ • Figure 3-2 
is a graph of the average number of aircraft in the interval [x, x -f- Ax) vs. 
x. It is assumed that the velocity is constant. For the example shown, leg 
1 and leg 3 are the most commonly entered legs. Note that the model 
correctly exhibits increased congestion in the traffic pattern as the touch- 
down point is approached. 


3.3 Nonhomogeneous Poisson Airspace Model 

The results of Sections 3.1 and 3.2 can be generalized by allowing the 
number of aircraft landings in the interval [o,t) to be a nonhomogeneous 
Poisson process with rate parameter 


X 


t 



v(a) da 


(3-20) 


where v(a) > o is a variable intensity parameter that accounts for dally 
fluctuations in air traffic density. That is. 


,n 


P[N(t) - n] - 



(3-21) 


Proceeding as before, assign probability to the arrival of an 

aircraft at leg m and denote by N^Ctjx) the number of aircraft entering the 

landing pattern at leg m in the t-second interval [t, t + t) . Thus, 

for the distribution of N (t;t), 

m 


P[N (t;r) = n ] = 
1 ra ’ m 



PlN m (t;T) 


n /N(t + t + t) - N(t + t) 
mm m 


n] 


P[N(t + t + r) - N(t + t) = n] . 
m m 


(3-22) 
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Average No. of aircraft in 
[x, x + Ax) at time t 
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This uses the law of total probability and exploits the fact that 

every aircraft entering the landing pattern at leg m must land t seconds 

later. Noting that N(t + t + t) - N(t + t) is Poisson with rate parameter 

mm 

+ t) “ + T - X £ + , eq. (3-22) can be evaluated as follows 

m m 


P[N (t;r) 
m 


n ] 
m 


m , \ n 

e CO 

n“n m v 

m 


U " p„> 

m 


n-n 3° ( t ; t + t) 

, m m 


-6(t;t m + t) 


nt 


n 


[P 


m 


(t»t m + t) ] 

n 1 

m 


m 


- p » 6<tit m + T) 


(3-23) 


Thus N m (t;t) is Poisson with rate parameter p m 3(t;t m + t) : 


\(t;x) : P[p m fl<t ; t m + r)] 


t + t + T 

/ m 

t + T 
m 


v(a) da 


(3-24) 


Proceeding, as in Section 3.1, it can be shown that N^(t;i), N^ftjT), 
■ Nj yI (t;t) are independent processes. Note here only that the rate 
parameters (i.e., average number of aircraft entries) for N^(t), 

N^(t) give the correct rate parameter for N(t). The property is pre- 
served that the average number of aircraft landings in the interval 

t + t^) equals the average number of aircraft entering the lauding 
pattern at leg 1 in the interval [o,t) plus the average number of aircraft 
entering at leg 2 in the interval t + t^ - t^) plus the 

average number entering leg M on the interval [t n - t w , t + t, - t w ) . To 

2 M 2 M 

show this 
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v(ot) da - + p 2 0(t;t 2 + - t 2 ) + 


i * * 


t + t. 


/ 


+ p m^ t;t M + fc l 


V 


M 




j-l 




(3-25) 


and therefore the averages are preserved as required. 

One can proceed as in Section 3.2 to determine the number of aircraft 
in the spatial interval [x,x + Ax) in a nonhomogeneous airspace. The 
results are direct extensions of eqs. (3-18) and (3-19). 

3.4 Multidimensional Extensions to the Basic Airspace Model 

Multidimensional extensions to the basic airspace model are obtained 
by considering processes of the form N(t) with t a multidimensional t-set. 
Then for example, one can define t = [o,t) X [o»x) X [o,y) X [ 0 , 2 ) to be a 
cross-product set consisting of temporal and spatial variables. The process 
N(t) becomes the number of aircraft contained in the region of space [o,x) X 
[o,y) x [o,z) during the time interval [o,t). If N(t) is a Poisson process, 
then the relevant probability statement is 

-Xt 

P[N(t) = n] = ~ e t (3-26) 

where X is now the average number of aircraft in the t-set t. The rate 
t _ 

parameter is related to a multidimensional intensity paramter, v(t). 
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describing the temporal-spatial Intensity of aircraft in the space: 


- 


/ 


da v(a) 


( 3 - 27 ) 


The intensity parameter v(t) is, in turn, chosen to describe the geometric 
properties of a typical aircraft trajectory in a pre-speclfied landing 
pattern. 
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CHAPTER 4 


MODELS FOR HAZARD PROBABILITY 

In this chapter a hazard Is modelled as the occupancy by two or more 
aircraft of a common region of airspace. The hazard model is tied to the 
basic airspace model of Chapter 3 and analytical expressions are derived 
for hazard probability. When aircraft trajectories are fixed-velocity and 
deterministic, the hazard probability can be readily evaluated. Numerical 
results Indicate that th parameter > with v the airspace intensity 

parameter and a/v 0 the so-called dwell time of an aircraft in a hazard-free 
interval, is a reasonable figure of merit for uncontrolled airspaces. Hazard 
probability is reduced by reducing £■ 

Some discussion is given to a hazard model that is applicable when 
aircraft trajectories are random. The resulting expressions for hazard 
probability can be evaluated for arbitrary random trajectories. 

4.1 An Occupancy Model for Deterministic Trajectories 

Assume all aircraft in the landing pattern of Figure 3-1 are flying 
deterministic, fixed velocity (v ) trajectories. In Section 3.1 it was 
shown that Q(x;t), the number of aircraft occupying the spatial interval 
[x,x + Ax), d <x<x+Ax<d . - , Is Poisson distributed: 

Q(x;t) : P[p] 


b - 


yAx 

v 

o 


E 


j-i 



(4—1) 


In this model, U is the average number of aircraft in the region [x,x + Ax) 
at time t. 

Now assume that the number of aircraft in [x, x + Ax) is 

Q(x;t) = q. It follows from the basic properties of the Poisson process 

that the unordered occupancy points for the q aircraft are uniformly 

distributed in [o,t). The ordered occupancy points x x. < x„ < . .. x < 

tH -L & 4 

x + Ax, where x^ is the location of the i aircraft, are distributed as 

the order statistics U^, . .., U ; that Is, 

q 27 
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f-(u) - 
U 


5i_ 

Ax q 


, x < u. < u„ < . . . < u < x + Ax 

f — I 9 n — 


0 , otherwise 


(4-2) 


where fjj(u) denotes the density function for U ■ U 2 U^) . 

A rather classical result that will be needed for our analysis of 
hazard probability is the probability that all inter-aircraft distances 
u 1+1 ~ u^ exceed a hazard-free distance a. This is the probability that 
x <_ u 1 < x + Ax - (q-l)a, u 1 + a < u 2 < x + Ax - (q-2)a, . . . ^ + a < 

u q < x + Ax. Denoting by P[H/Q(x;t) » q] the probability that the q-1 
inter-aircraft distances all exceed the hazard distance a (the choice of 
notation will be clear shortly) , then 


P[H/Q(x ; t) 



x+Ax-(q-l)a 

/ du i 
x 


x+Ax-(q-2)cx 

/ du 2 ' ' 
u l + “ 


I 4 


x + Ax 



(4-3) 


Letting “ (u^, - x)/Ax, 


P[H/Q(x;t) = q ] = q! 


l-(q-l)a/ Ax 
/ dY i 


o 


l-(q-2)ct/Ax 1 

/ dY 2 ■••/ ' 

Yj+a/4x Vl 4 ”* 


(4-4) 


Tills can be simplified (ref. 10) yielding 


P[H/Q(x;t) -q] 


1 , q _< 1 

[i - (q-l) ]“ , 1 < 4 <[f +l] 


0 


otherwise 


(4-5) 


2-8 



where [Ax/a + 1] denotes the largest integer less than or equal to Ax/a + 1. 
The parameter Ax/a is the number of contiguous hazard-free intervals in the 
interval [x, x + Ax). 

This yields a model for hazard probability for uncontrolled airspaces 
of the form depicted in Figure 3-1. A hazard is defined to exist in 
the airspace interval [x, x + Ax) at time t if two or more aircraft occupy a 
common region of airspace. If the common region of airspace is taken to be 
an interval of length a, then the result of eq. (4-5) can be used to model 
the probability that no hazard (H) exists in the airspace [x, x + Ax) at time 
t, given that q aircraft occupy the airspace Interval [x, x + Ax). The law 
of total probability then yields the following model for the probability that 
a hazard (H) exists in the airspace region [x, x + Ax) at time t: . 

P [H] « 1 - P[H] 


“ 1 - P[H/Q(x;t) = q] P[Q(x;t) = q] 

q=i 


r^+ 1 ] 

= 1 - e _1J [1+y] - e ~ W iT (1 " (q_1) Ax )q <4_6) 

q=2 


Equation (4-6) is a basic expression for the probability of a hazard. In Fig- 
ures 4-1 and 4-2, P[H] is plotted versus the rate parameter p and parametrized by 
Ax/a, the n' iber of contiguous, r.on-overlapping, hazard-free regions in 
[x, x + Ax), The parameter [Ax/a 4- 1] is the maximum number of aircraft 
that could, in principle, be placed on the Interval [x, x + Ax) without 
creating a hazard. The results of Figure 4-1 indicate that the hazard 
probability is small when p, the average number of aircraft to be found in the 
airspace interval [x, x + Ax) , is very much smaller than the number 
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Figure 4-1. Hazard Probability vs. Average Number of Aircraft. 
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Figure 4-3. Random Aircraft Trajectories. 
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of contiguous, non-overlapping, hazard-free regions in [x, x + Ax). That is, 
for small PfH), 


p « Ax/a . (4-7) 

Equivalently, a figure of merit can be defined as 


5 ■ P (Ax/a) ^ 



requiring f; <ai 1 for law hazard probability. The parameter «/v q is the 
time it takes a constant velocity aircraft to travel one hazard-free distance. 
Thus, for £ Pj = 1 (worst-case), the figure of merit £ is the average number 
of aircraft that enter the airspace in the time it takes an aircraft to 
traverse one hazird-free distance. From Figure 4-2 it is found, for example, 
that £ on the order of £ “ 0.01 ensures P[H] on the order of P[H] * 0.01. 

4.2 An Occupancy Model for Random Trajectories 

When aircraft trajectories are random rather than deterministic, 

Q(x;t), the number of aircraft occupying the airspace region [x, x + Ax) at 
time t, is no longer Poisson distributed* and the previously derived results 
for hazard probability must be modified. To achieve this modification, 
consider the diagram of Figure 4-3. Several nominal, fixed-velocity, 
deterministic trajectories are Illustrated, along with random trajectories 
that deviate from these nominals. The trajectories terminate at the 
Poisson event times t^, X 2 » ..., T n corresponding to the homogeneous process 
N(t) that characterizes the number of aircraft landings in th^ time interval 
[o,t). Note in the figure that the first aircraft to enter the airspace 
enters at leg 2, whereas the first aircraft to land enters at leg m. 
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It would be desirable to evaluate the probability that, given Q(x}t) - q, 
no two aircraft occupy the same region of airspace. Then the procedure of 
Section 4.1 could be used to evaluate a hazard probability. The difficulty is 
that even when the nominal, deterministic trajectories are separated by a 
distance a, the corresponding random trajectories may come well within a 
of each other. This effect is illustrated in Figure 4-3 where it is shown 
that two random trajectories in the region [x, x + Ax) at time t come 
within a (and therefore constitute a hazard) even though all nominal 
trajectories are separated by more than a. 

These difficulties imposed by the random trajectories can be overcome by 
requiring the nominal trajectories to be separated by a plus a random slack 
variable that guarantees the corresponding random trajectories are also 


separated by a. This strategy is illustrated in Figure 4-4. Trajectories 
j and j-1 (with landing timeB x^ and x^ ^) are required to be separated by 


the distance a + h. where 

j 


h j ’ ' Vi 


(4-9) 


and 3^ is the random perturbation of the j trajectory from its nominal 
trajectory at time t. 

Now consider the ordered occupancy points x £ x^ < x 2 < • • • £ x + Ax 

for the q nominal-trajectory aircraft located in [x, x + Ax), These occupancy 

points are distributed as the order statistics U. , U_, ..., U , as shown in 

1 / q 

eq. (4-3). The probability that the j-th and (j-l)-st aircraft are separated 
by a + hj for j = 2, 3, . . . , q is the probability that u_j - u . ^ > a + h^ . 
Equivalently, this is the probability that all random trajectories are 
separated by at least a. For given h^ , h,,, h 3 , ..., h^, the appropriate 
probability statement is 


35 



P[H/Q(x;t) 


x + ix 



X 


(q-l)o 



(4-10) 
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x + Ax - (q-2)a - ^ h. 

j*3 -> 


Ax + x 


du, 


2 ‘ 


+ a + h 2 


/ 


du 


u , f a + h 
q-1 q 


Letting Yj 


(uj - x) / Ax , 


1 - 


P[H/rj(x;t) = 


q» h 1# ... h ] = qj 


/ 


q 

(q-!) Ax + Ax ^ h j 


dy. 


1 - <«- 2 > 5 


/ 


- £ 

j=3 J 


. a_ h 2 
Y 1 Ax Ax 


dy 2 ••• j 


f 


(4-11) 


, dy 

y *i- + !a n 

f q-l ' Ax Ax 


After simplification 
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?[H/Q(x;t) - q, h v .... 


I I , a 1 

[1- <q-D fj- (6, -B 1 )] q , Kq <[f +l] • 

0 , otherwise (4-12) 


Defining the random variable y = 
is independent of q, for 1 < q <_ 


B - and assuming the perturbation B 

^A-v ^ 


[ — + 1 ], 

O 


P[H/Q(x;t) 



“ PfH/Q(xjt) 


q> y] 


tl - (q-« fc - Tl" • 


(4-13) 


The result here is very similar to the result of eq. (4-5). However, to obtain 
the unconditional probability of a hazard, P[H], in this model, averaging 
over the joint distribution of Q(x;t) and h is necessary. Denoting the 
distribution of y by F , 


P[H] = 1 - P[H] 


= 1 - 


CO 


E 

q=l 



P[H/Q(x;t) 


q, y] P [Q - q] 


(4-14) 


- 1 - e 


V 





This is the basic result for hazard probability when aircraft trajectories 
are random. Note that by using the binomial expansion • for (a + b)^, 
eq. (4-14) can be written in terms of the moments of yi 
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P[H] 


[1 + M] - 


\& + ll 

La J 

2-i 1 

q*2 


qi 


l 

Ax q 


g (2) tAx - (q -W«] k / 


dV y (-Y) q " k 


« 1 - e" 1-1 [1 + y] 




[Ax - (q-l)a] k 
(q-k) I kl 


(-l) q ~ k in. 


q-k 


where (4-15) 

Vic - ElY q “ k ] 


Thus, by Bimply knowing the first [— ■+ lj moments of y» one can in principle 

compute the hazard probability. When the slack variable y is normally 

distributed, say with mean zero and variance cr^, then m . » 1*3 ... (q-k-l)a^ q k ^ 

q-K 

for q-k even and zero otherwise. Then eq. (4-15) can be easily evaluated. 

4.3 Multidimensional Extensions 

Multidimensional extensions of the results given in Sections 4.1 and 4.2 
are, in principle, straightforward. One simply considers multidimensional 
Poisson processes N(t), with t the multidimensional t-set discussed in 
Section 3.3 and examines the probability that aircraft are separated by a 
suitable function of the multidimensional variable t. In this way geometrical 
considerations for hazard regions can be Included and more complicated 
trajectories analyzed. The results for hazard probability will be multi- 
dimensional analogs of eqs . (4-6) and (4-14) . 
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CHAPTER 5 


SIMULATED AIR SPACES 
5 . 1 General 

A preliminary study has been conducted on data previously collected 
by NASA Wallops personnel at three uncontrolled airports in Maryland. The 
data consist of single radar tracks of general aviation aircraft in an 
uncontrolled environment. Consideration is given to ways in which these 
data can be used to evaluate alternative patterns and procedures in such 
environments. To evaluate collision risk in an air traffic environment 
one approach is to analyze an environment which provides representative 
space-time relativ'nships between aircraft for given patterns and procedures. 
The real data do provide spatial relationships which are representative of 
pilot adherence to the current procedure for runway approach. 

To evaluate current procedure one obvious method is to assign some 
time correlation to the real data tracks and then sample this finite set of 
position-time situations. Sampling could be done on a "with replacement” 
basis to provide an indefinitely large (but finite) set of "situations." 

The method is direct and should not be difficult to implement but it pre- 
cludes tracks other th*m the specific ones in the data base (tracks which 
are plausible but did not occur) from being analyzed. This could conceivably 
introduce a bias in the results. Furthermore, since it is desired to develop 
the capability to evaluate other patterns and procedures besides those which 
gave rise to the tracks in the data base, a simulation model to generate 
tracks for analysis is desirable. Another consideration involves the fact 
that the real data base is not presently available in a form which Is 
suitable for computer analysis. Thus to even validate, for example, the 
statistical occupancy model evaluation discussed in Chapters 3 and 4, it 
is necessary to have some simulation model to provide track data. 

The simulation model to be discussed here is a generating approach in 
the sense that tracks are generated which are plausible in light of the 
data, but which are not themselves in the data base. This approach would 
regard those tracks in the data base as a sample from which population 
distributions would be derived. The tracks to be used in collision risk 


39 



analysis would then be generated using some Monte Carlo sampling from such 
distributions. This approach minimizes the danger of bias, is readily 
adaptable to different patterns and procedures, and can be implemented 
very easily once a method is defined to characterize air traffic. 

One of the first issues addressed in this study to provide a basis for 
a simulation model was how traffic flows can be characterized. This issue 
Is important in identifying different traffic flows and in the validation 
of computer models. The resulting approach presented here as a recommendation 
of this study is based on the major segments of tracks flown, classified by 
type of aircraft, and each aircraft entry point into a defined airspace. 

Each segment is characterized by a few basic flight variables such as headings, 
segment lengths, turn radii, speed, altitude, etc. This approach seems 
feasible from the standpoint of the capabilities of the data base and appears 
adequate for the purposes of the analysis considering precision requirements 
and cost of analysis. 

Sections to follow provide a description of a prototype digital simulation 
model which has been used to generate traffic for the existing environment. 

A capability is described to model simultaneous tracks in a multiple aircraft 
environment designed to represent an uncontrolled terminal area airspace. 

The model should be applicable in an existing situation as well as in a 
hypothetical environment which can be reasonably postulated. Central to the 
development of this capability is a means of generating a virtually 
unlimited number of tracks in the environment to be modelled. The descriptions 
are not complete but are intended to convey an idea of how the simulation 
model would work and how it can be used to analytically evaluate the traffic 
patterns from a mid-air collision hazard standpoint. 

Data requirements fur the prototype simulator are discussed. 

Essentially these are empirical distributions of path lengths, headings, and 
turn-radii of legs that make up a track. These requirements can be readily 
met by the data retrieval and analysis programs now being developed by NASA. 
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5.2 Characterization of Air Traffic Flows 


For purposes of this development the uncontrolled airport airspace 
will not be the entire "Airport Traffic Area" conventionally defined as 
that airspace within a 5 statute mile (8 km) radius of the center of the 
airport up to but not including 2000 ft (609.6 m) altitude. The available 
track, data in the NASA Wallops data base drops off rapidly beyond a 3 n.mi. 
radius (5.6 km) of the radar site. Furthermore the main feature of 
Interest is the prescribed landing approach pattern which is usually flown 
within this 3 mile radius of the runway. Therefore, a perimeter of 3 n.mi 
(5.6 km) radius from the runway threshold will be defined to encompass the 
airport airspace in the pursuing analytical discussion. A coordinate 
system with origin at the threshold and the y-axis coincident with the 
runway centerline as shown in Figure 5-1 will be used. This choice of 
coordinate system should not prejudice the main concepts and results to be 
presented. 

The digital traffic generation model described represents a simple 
yet plausible characterization of the uncontrolled environment. Aircraft 
movements are described in terms of a few basic flight variables. Tracks 
are defined in terms of major component segments and the general environment 
is defined in terms of traffic types and the prevailing patterns and 
procedures. 

5.2.1 Relevant parameters In simulated traffic studies entry of 
aircraft into the defined airspace is often the fundamental "driving event," 
Hence, a logical first variable of characterization is the time between 
entries.* This includes arrivals to the perimeter of aircraft coming in 
to land, or flying through, as well as aircraft presenting themselves at 
runways for take-off. Time between entries is a variable that often has 
diurnal and other longer term variations which can be accounted for as 
necessary. 


*Vlrtually all of the variables described are best thought of as random 
variables, and when considering a variable for characterization its 
empirical frequency distribution should be used. The totality of such 
distributions will be said to characterize the traffic. 
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Next, a classification of various types of driving events and the 
frequency with which they occur should be considered. The simplest classi- 
fication, which should be adequate here, is "arrivals,” "departures" and 
"fly-bys." Arrivals and departures are commonly the greater portion of 
airport traffic relative to fly-bys and in the long run must be of equal 
frequency for most analytical purposes. Thus it should be enough to state 
the fly-bys as a fraction of arrivals, or perhaps airport operations. 

Having characterized traffic by its inter-event time and by type of 
event, further characterization is bast done separately for the three types 
of events. In what follows, only the characterization of arrivals is 
considered. This operation is somewhat more complex than the others and 
is perhaps of greater interest for collision hazard studies since propor- 
tionally, arrivals spend more time in the airspace than other types. In 
any case, the basic approach to traffic characterization would be similar 
for the other types of operations also. 

The arrival traffic flow within a given airspace depends greatly on 
the type of aircraft, the point from which it enters and the runway on 
which it lands. Therefore, the distribution of arrivals by type along 
the perimeter and the distribution of runway usage are important items 
for characterization of traffic. 

For a given type of aircraft, entry point, and runway, the flow 
patterns of traffic are governed by a variety of factors such as rules 
and procedures, weather conditions, presence of other aircraft, individual 
characteristics of aircraft and pilots, etc. No attempt will be made to 
characterize traffic flows by any specific factor in the above group. It 
will be enough to accept empirically whatever flows do evolve under the 
influence of all of these factors, and to cry to characterize flows on 
the basis of major segments in the flight path. 

In order to define a segment, it is necessary to recognize six modes 
of flight: straight climb, straight descent, straight-and-level, climbing 

turn, descending turn^and flat-turn. Then a segment is defined as the 
portion of flight path over which a particular mode is maintained; the start 
of a new mode marks the start of a new segment. This essentially simple 
concept poses two practical problems. One is the difficulty of identifying 
a sharp changeover point from one mode or segment to another; the other is 
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the distinction between the "fundamental" mode and the "random" deviations 
from it. These problems will not be completely and quantitatively resolved 
here. The model will be designed to recognize no more than (k + 4) segments 
to any path, where k is a number that depends on the entry point of the 
arrival and represents some reasonable number of segments necessary to fly 
the recommended pattern for the airport. For example, in Figure 5-1, for 
arrivals from the sector shown, k may be taken to be 7: Inbound, turn on 

downwind, downwind, turn on base, base, turn on final, and final. Hence, 
the assumption can be made that no path through this sector need have more 
than 11 segments. If more than 11 segments seem to exist in fact, then it 
will be assumed that some of the apparent segments are indeed "ripples" and 
will be combined into a total of 11 or fewer segments. 

The above scheme is very helpful in characterizing traffic flows because 
it provides a means of classifying flight paths — by the number and sequence 
of segment types. The number of different types of flight paths (for a given 
aircraft type from a given sector) is likely to be quite manageable, because 
many segment types are in practice incompatible with one another, and it is 
possible to specify the frequency with which different flight path types 
are found at a given airport. Thus, the characterization of traffic at an 
airport should include a specification of empirical distribution of flight 
path types. 

Finally, for each type of flight path, the segments need to be char- 
acterized by distributions of a few basic variables: headings, speeds, 

altitudes, decelerations, descent rates, and turn radii. 

The above scheme represents a convenient way of characterizing traffic 
at an airport. In fact, it leads directly to the method of track regenera- 
tion and traffic simulation to be discussed in following subsections. 

5.2,2 Track generation .— In the generating approach, the track is 
developed point by point and segment by segment. The first step is to 
locate the arrival on the perimeter. This is done by sampling from dis- 
tributions for entry point, initial speed, initial altitude, and initial 
heading for the particular segment and type of aircraft involved. It may 



be that entry speed, altitude and heading will be dependent upon entry 
point and this dependence would need to be quantised before defining 
distributions. This is markedly true for Initial headings, as m:.y be 
seen from the data base, and the dependence can be identified from a 
plot of initial headings versus entry points. 

Having established the arrival on the perimeter, the next step is to 
generate what will be called the initial leg; this is the segment flown 
immediately after entering the perimeter. Again, the nature of this 
segment (with respect to the 6 modes mentioned earlier) depends very much 
on the entry point. From a brief observation of the data it seems there 
are some entry points with a large proportion of straight-and-level initial 
legs, while other entry points exhibit a large proportion of low-bank-angle 
descending turns directly onto the final course. The proper distribution 
for the particular entry point needs to be set up and sampled. 

Assuming for the sake of further illustration that Lhe initial leg 
is determined to be straight and level, a distribution of leg lengths for 
the particular aircraft type and sector is sampled. Trom the Information 
developed thus far it is then possible to generate the path of the arrival 
in space and time from the perimeter up to some point within the airspace. 

Next, the nature of the second segment needs co be determined by 
sampling from an appropriate distribution. Assume it is a descending 
right turn. Then the rate of turn, descent altitude and final heading are 
determined by sampling, and the generated path is extended to include the 
turn. 

This turn may or may not put the aircraft in the pattern; it depends 
on the entry point. Assume for the sake of discussion that it puts the 
aircraft on downwind. Then the next step would be to sample proper dis- 
tributions for downwind headings and other characteristics of the downwind 
leg, so that the path can be generated up to the point where the base turn 
begins. 
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The number of segments can vary from 1 (for sweeping turns directly 
onto final) to about 12 (for aircraft beginning with upwind legs and 
executing a full pattern) depending on point of entry; similarly, the 
nature of the segments will vary, depending 0 n type of aircraft. All 
these relationships can be gleaned from the data and programmed into the 
logic of the track generating routines. The feasibility of this approach 
has been checked out by some preliminary test programs, and the approach 
seems quite capable of generating tracks which are a faithful reproduction 
of reality- — the quality of the reproduction being determined by the 
level of detail employed in the set of sampling distributions used. 

To evaluate collision risk in an air traffic environment it will be 
necessary to create an environment which provides some time-space relation- 
ship between aircraft which can be representative of an actual environment. 
The simulation approach as described thus far provides a means of character- 
izing flight paths such that a path generator algorithm can recreate single 
track data. 

For the model to adequately represent the air traffic environment it 
is necessary to provide time correlation between aircraft cracks. The 
method proposed and used to illustrate how the model can provide baseline 
evaluation of collision risk involves assignment of entry times to single 
tracks generated using Monte Carlo sampling. The same procedure could be 
used to reconstruct a meaningful time correlation between individual tracks 
in the NA? \-Wallops data. Thu? each track generated in the airspace will 
have an entry time assigned (randomly) which then defines a space-time 
history to that track with an assumed velocity for the aircraft type 
corresponding to the track. 

5.2.3 Model validation . - Once the model is implemented there is a 
need to provide some means of validation. One method is actually a two- 
phase approach. A first validation can be done on all of the same variables 
that are used to characterize the traffic. Distributions would be compiled 
on these variables from a large number of generated tracks, and compared 
against corresponding distributions from data base tracks of the same air- 
craft type and entry point. If the distributions match to some stipulated 
goodness of fit, the generated tracks may be considered valid. Because of 


46 



the close relationships between the distributions generated, the distribu- 
tions sampled from and the distributions developed from the data, the above 
test by Itself does not provide adequate reassurance as to the validity of 
the generated tracks. Hence, a second-level validation is proposed which 
does not involve distributions used in track generation. These could bo 
cross-sectional distributions of the type investigated in previous NASA 
studies. A number of planes can be selected and cross-sectional distribu- 
tions developed from the data base as well as the generated tracks. If the 
distributions match to some stated degree of r, .t, the generated tracks may 
be considered further validated; if not, more detail needs to be incorporated 
in the set of generating distributions. Of course, all this needs to be 
done with due regard to sample sizes, sampling errors, and all the other 
requirements of statistical testing. 

The following sections describe the path generation simulation model 
as applied to two of twelve possible entry sectors. The algorithm for 
assignment of interarrival times is included in this discussion. In Chapter 
6 model validation is discussed and illustrated. 

5.3 Air Traffic Simulator for the Uncontrolled Environment 

In Section 5.2 air traffic path characteristics ware discussed in 
terms of variables which can provide for computer generation of aircraft 
tracks. The procedure for generating tracks was described in general. Ap- 
pendix A contains a detailed description of how entries from two of the twelve 
possible entry sectors are assigned paths within the airport airspace. 

Path segments are generated by random sampling from parametrized distri- 
butions. These parameters are defined so that the distributions can be 
changed to conform with distributions that could be obtained from real 
data or might be analytically determined from chosen pattern characteristics 
and procedures prescribed for a given uncontrolled airport. 

Once this ability to generate tracks is developed as described 
in Appendix A it is possible to design a computer simulation model to study 
the collision hazard in a given environment. This section describes the 
model and how it operates. 
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Tha computer modal has the following characteristics: 

(1) It assigns an arrival time for each aircraft (if time is to be a 
factor in the simulation). Tha aircraft arrival time is based on the as- 
sumption that the inter-arrival times are exponential. The following equation 
was employed to generate exponential inter-arrival times: 

TAU * - j £n(u) 

where A ■ number of aircraf t/unit time, 

u = uniform random number. 

(2) The flight paths are generated as described in Appendix A. Points 
along the flight path are calculated based on a specific time interval and 
velocity. Each flight may be defined in both time and space. 

(3) The amount of computer storage required to maintain the flight data 
was kept to a minimum by examining flight paths in sequential pairs. To 
illustrate, the first flight pair would be AC1 and AC2 where AC1 was the 
prior flight and AC2 is the current flight. Once data are obtained from the 
two flights and the statistics updated, another flight is generated and AC2 
becomes the prior and AC3 becomes the currert. This procedure simplifies 
the calculation of the time parameter in that the prior flight will alwa; a 

be entering the airspace at time zero and the current flight will be entering 
at time TAU. 

The computer model consists of five logical units as shown in Figure 5-2. 
The Initialization unit handles variable assignment and sets specific para- 
meters by user input. It Initiates the entire air pattern simulation. The 
Sector Selector handles the selection of the appropriate sector from which 
an aircraft will be entering the airspace. It also tallies the total number 
of flight paths simulated. The Path Generators are sets of programs (one 
for each sector) that operate flight paths based or. the sector entry tra- 
jectories (see Appendix A) . The Statistics Generator maintains desired 
histograms and updates the recursive mean and sums of squares functions. The 
Output unit handles the output of the desired histograms and related means 
and standard deviations. 
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There are currently five main programs that comprise the Air Pattern 
simulator. These are outlined in Appendix B along with the logic unit 
association. Individual flowcharts for each program are also given in 
Appendix U with a brief description. 

For collision risk analysis the simulation model output consisting of 
prescribed histograms and statistics is used as input to special analysis 
programs. The types of output available and some of the analysis methods 
used are discussed in Chapter 6 along with presentation of selected results. 
The types of output available include that necessary to provide validation 
of the model. In the computation and interpretation of hazard statistics 
proper account must be taken of the fact that the tracks generated do not 
reflect any avoidance maneuvers. This is a basic limitation of the simulator 
which is really imposed by limitations in the real data. The data do not 
provide information concerning what parts, if any, are the results of pilot 
maneuvering to avoid a potential conflict situation. This limitation could 
be presumably overcome by Including decision making algorithms in the 
simulation model. In doing this ^t would be desirable to obtain surveil- 
lance type data which could give a composite picture of the traffic at all 
times and the maneuvers performed in conflict situations. 
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CHAPTER 6 


STATISTICAL ANALYSIS OF SIMULATED AIRSPACES 

6*1 Histograms of Spatial Trajectory Deviations 

6.1.1 Numerical results .— In order to verify that the trajectory model 
could generate valid flight profiles (in at least two dimensions), 500 
aircraft were cycled through the approach pattern and statistics were 
calculated at downwind 1 (DWl) , downwind 2 (DW2) , and base plane (BASE) 
locations (see ref. 2), In addition to histograms of distance deviation, 
histograms of distance difference and time difference of sequential air- 
craft pairs were generated. All data presented in this section were 
generated with the arrival rate parameter selected at 100 aircraft per hour. 
The results may be generally summarized as follows: the distance histograms 

appear nearly normal with slight skew at the longer distances, the distance 
difference histograms appear to resemble a one-sided normal curve, and the 
time difference histograms appear exponential as would be expected. These 
data are included in Figures 6-1 through 6-9. 

6,1.2. Comparison with published data .— Figure 6-10 shows selected distance 
deviation da^a from actual aircraft observations (ref. 2). Only distance 
data at SW1, DW2 and BASE planes are shown. Notice that in general the 
histograms are similar to those generated by the trajectory model with 
the exception of the increased skew at larger deviations. This difference 
in skew may be attributed in part to the random number generator in the 
trajectory model and in part to the skewness introduced by single engine 
high wing aircraft in the observations (see Figure 6-10d) . Table 6-1 shows 
a summary comparison of first and second order statistics between actual 
and simulated distance deviations. Agreement is good at the DWl and DW2 
planes but tends to appear suspect at the BASE plane. This can perhaps be 
attributed to trajectory model constraints in the base leg turn and the fact 
that approaches directly into the base leg arc not roneidered in the current 
simulation. 
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Figure 6-1. Distance Deviation - Downwind 1 Plane. 
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Figure 6-2. 


Distance Difference for Sequential Aircraft, 
Downwind 1 Plane. 
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Figure 6-3. Time Difference for Sequential Aircraft 
Downwind 1 Plane. 
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Figure 6-4. Distance Deviation - Downwind 2 Plane. 
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Figure 6-6. Time Difference for Sequential Aircraft, Downwind 
2 Plane. 
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Figure 6-8. Distance Difference for Sequential Aircraft, 
Base Leg Plane. 
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Figure 6-9. Time Difference for Sequential Aircraft, Base 
Leg Plane. 
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Figure 6-10. Published Data (ref. 2) from Actual Aircraft Observations 
in Uncontrolled Terminal Area. 







Table 6-1. Comparison of Actual and Simulated Statistics 
at DW1, DW2, and BASE Planes. 



Actual (ref 

.2) 

Simulated 



mean 

std dev 

mean 

std dev 

DW1 

-5440' 

2402' 

-5200' 

1162’ 

DW2 

-5600' 

780' 

-4869' 

1561' 

BASE 

-5535 

552' 

-3414 

1986' 
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It Is apparent from the preceding discussion that the current simulation 
can to a limited degree describe an uncontrolled terminal area in terms of the 
types of data previously published. It Is felt that additional flexibility, 
when Incorporated, will allow the trajectory model to be utilized in concert 
with the theoretical model discussed in previous chapters to assess hazard 
probabilities . 


6.2 Estimated Hazard Probabilities 

The trajectory model can be used to generate an estimate of hazard 
probabilities by observing the simulated airspace on an jircraft pairwise 
basis. If the simulation Is exercised in such a manner that the separation 
distance of two aircraft occupying the area jointly is indexed versus time, 
and if this is performed over all aircraft pairs, then several outputs 
be. ome available. 

The least complex of these Is a histogram of separation distance 
over the entirp ensemble of aircraft pairs, While being the simplest output, 
it perhaps is also the least meaningful in predicting or assessing a hazard 
except to give an intuitive feel for the value of distance to closest 
approach which would (or should) be used to estimate hazard status from the 
other available outputs. 

Figures 6-11 through 6-13 show histograms of separation distance for 
intensity ^ 'ameters of 50, 100 and 200 aircraf t/hour. The shape of the 
histograms for 50 and 100 aircraft/hour are similar in that a predominance 
is indicated about line number 20 (3000 ft) and that a relatively uniform 
distribution of distances is observed elsewhere. The histogram for 200 
aircraft/hour begins to take on fairly uniform characteristics over its 
entire range and apparently indicates that a saturation situation exists 
and intensity parameters in this range should not be expected to generate 
realistic results. It would appear that minimum separation distances of one 
to three thousand feet are proper parameters for evaluating hazard probabilities 
(within the influence imposed by the limited capability simulation) . It should 
also be remarked that distances greater than ten thousand feet were deleted 
from consideration. 
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Figure 6-11. Distance Histogram for 50 Aircraf t/liour 
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A more meaningful assessment of the terminal area environment is what 
can be defined as "pilot workload." This is manifested as a histogram of 
the total time two aircraft spend at a separation distance less than some 
prespecified distance. This represents that distance at which a pilot would 
be concerned with the decision to make a see-and-avoid maneuver. 

Figures 6-14 and 6-15 show histograms of encounter time length for an 
intensity parameter of 100 aircraft/hour and for minimum separation distances 
of 2000 feet and 1000 feet. Notice that the histograms appear exponential 
as would be expected from both the nature of the random variable and the 
density function assigned to inter-arrival times in the simulation. In going 
from a separation distance of 2000 feet to one of 1000 feet, the number of 
entries in the histogram decreased from 114 to 55. The histogram for 100. ft 
is thus observed to be more erratic than the one for 2000 ft. 

An alternate approach to representing "pilot workload" is to compute 
the cumulative time in an encounter normalized to the total time spent in 
the terminal area. Figures 6-16 and 6-17 show histograms demonstrating this 
type of data. Notice that as the separation distance parameter is decreased 
to 1000 feet, the prominence near bin 20 moves toward the origin. This is 
because, as the separation parameter is decreased, the number of aircraft 
which have the opportunity to spend a long time in an encounter is reduced 
(for a fixed number of aircraft) and more short encounters are produced 
(i.e. for zero leet separation parameter, there would be zero encounters 
while for infinite separation parameter, all aircraft would be in an encounter 
100% of the time). Notice that for 2000 feet separation, 2% of the aircraft 
are in an encounter 100% of the time in the terminal area (see bin 51). 

Table 6-2 shows a summary of the histogram data including means and 
standard deviations. Notice that as the Intensity parameter varies for a 
given separation distance, the number of the entries for each histogram varies 
accordingly but that the mean and standard deviation do not vary appreciably. 
This is likely due to the fact that the 3ame paths are generated for each 
case and suggests that 100 aircraft are sufficient to generate good sample 
means and standard deviations. Notice that as the separation distance 
parameter varies, the means and standard deviations vary accordingly. 
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Histogram of Percentage of Time in an Encounter, 
100 Aircraft/Hour and 2000 Feet Separation. 
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Figure 6-17. 


Histogram of Percentage of Time in an Encounte 
100 Aircraf t/Hour and 1000 Feet Separation. 



Table 6-2. Histogram Summary 
Encounter Length Histograms 

(Percentage of Aircraft having a Given Length of Tima in an Encounter) 
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Another meaningful output available from the trajectory model is the 
number of Independent times, "encounters," two aircraft aro within a pre- 
specified distance. This is independent of the length of time spnnt in this 
situation 'ind likely becomes the output most closely related to hazard 
probability. For example, if one incorporated distance, rate, relative 
heading, velocity, and binoculars, this output could be used in a Monte 
Carlo fashion to produce estimated probability of mid-air collision in the 
same sense that it is done theoretically in previous chapters. 

To demonstrate the potential for using the simulation to examine hazard 
probabilities, the number of independent "encounters" was generated for minimum 
separation distances of 1000, 2000, and 3000 feet and for intensity parameters 
of 50, 100, and 200 aircraft per hour. The results are shorn in Figure 6-18. 
Referring to Figure 4-1 and invoking the rationale that: 

(1) AX be extended to the distance an aircraft travels from entering 
the airspace to touchdown, 

(2) a corresponds to the minimum separation distance, 

(3) the intensity parameter can be interpreted as a measure of 
average number of aircraft in the airspace, and 

(4) the number of encounters for for 100 aircraft, the average number 
per aircraft) can be interpreted as a measure of hazard probability. 

One can observe that the simulation does indeed produce data In consonance 
with theoretical predictions. 

Notice that the curves of Figure 6-18 tend to "flatten out" at high 
intensity parameters as opposed to the behavior observed with the large average 
number of aircraft in Figure 4-1. This is a direct result of the "saturation" 
mentioned previously as the intensity parameter is allowed to Increase. In 
fact, for very large values the curves would peak, then diminish, becoming 
asymptotic to 100 encounters. That is, as the number of aircraft increases, 
a point is reached where a give", aircraft is always in an encounter. Thus 
one obtains one and only one encounter per aircraft over 100 aircraft, 
producing 100 encounters independent of parameter value. This indicates the 
simulation is valid only over a practical range and is not useful in 
examining limiting conditions. 
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Figure 6-18. Number of Encounters vs Aircraft /Hour with 
Minimum Separation Distance as a Parameter. 
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CHAPTER 7 


A NON-CLASSICAL MODEL FOR RANDOM AIRCRAFT TRAJECTORIES 

The simulated trajectory model of Chapter 5 re present.* one reasonable 
approach to the simulation of random aircraft trajectories that enter a 
landing pattern under the Poisson regime. There are other models, such 
as the Gauss-Markov dynamical models popularized in the modern contrul 
literature, that one might consider. The application of such models to 
trajectory simulation and analysis is a well-developed, well-published 
topic. 

In this section departure is made from the Gauss-Markov dynamical model to 
consider Poisson-driven dynamical systems. The purpose is merely to 
suggest that there are models other than Gauss-Markov models that deserve 
consideration. The actual choice of a suitable model must be based on 
careful statistical analysis of competing models and comparison with 
measured data. 


7.1 Poisson-Driven Dynamical System 

The construction of a Poisson-driven dynamical model for random air- 
craft trajectories proceeds as follows. The random telegraph wave (ref. 11) 

T(t) = (-1) NC0 (7-1) 

with N(t) a Poisson counting process is defined as the control of the air- 
craft. This control is fixed at + l r between event times of the Poisson 
process as shown in Figure 7-1. The simplest of all Poisson-driven 
dynamical models for random trajectories is then obtained by letting y 
denote the iss-path deviation of a trajectory from the nominal landing 
pattern and modelling y as 

t 

y t = T(t) ; T(t) = J (-1) NC0) da . (7-2) 

o 
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Thus y Is an Integrated version of the random telegraph wave that represents 
fixed-velocity aircraft trajectories which deviate linearly from the nominal 
trajectory until they are corrected at a Poisson event time. A typical 
trajectory is illustrated in Figure 7-1. The process T(t) is often called 
a "random time" (ref, 11). 

A refinement to eq. (7-1) is obtained by introducing a damping term 0 
and writing y as 




- BU-tp 

e 



- 6Ct 
e 


a) 


(- 1 ) 


N (a) 


do 


(7-3) 


The corresponding random differential equation Is 


dy fc + 0y t dt = dT(t) 


(7-4) 


A typical trajectory for this model is also illustrated in Figure 7-1. The 
model can, of course, be generalized further to account for inertial effects. 

7.2 Density Evolution for Poisson-Driven Dynamical Systems 

The statistical characterization for y fc is complete when the joint 


density function for y , y.. , . . . , y is known fo ■ all finite t-sets 

C 1 u 2 ^n 

(t 1 , t_, . . . , t ) . In this section consider the determination of the partial 
X i n 

differential equation that characterizes the first-order density function, 
denoted f(y;t), for y t - In this derivation the increment dl(t) plays an 
important role. 
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Typical sample functions for dT(t) 


T(a) da are illustrated in 


t 4* dt 



Figure 7~2 for various event times t < a < t 4- dt. 
is 


The length of an increment 


dT(t) 


■/ 


t + dt 




dt 

■/ 


( -l) N(t + u) du 


(7-5) 


(dt - 2a)(-l) N(t) , o < a i dt 


Given that a Poisson event has occurred in the interval [t, t + dt) , the 
event time a is uniformly-distributed in [t, t + dt) . It can the 1 ' be shown 
that the distribution function for dT(t) is 


F dT <Y) = P[dT(t) < y] 

0 , y < - dt 

[1 - Xdt] P q 4- ~ (y + dt), -dt _< y < dt (7-6) 

3 , y > dt 

where P q is the probability that N(t) is odd and X is the intensity parameter 
for JJ(t). This distribution function, illustrated in Figure 7-3, has a jump 
of size [1 - Xdt] at y = -dt and a jump of size [1 - Xdt] [1 - P q ] at y s dt. 
The result simply says the Increment dT(t) has length -dt with probability 
[1 - Xdt] P q , which is the probability that N(t) is odd and no events occur 
in j.t, t 4- dt) , and length dt with probability [1 - *dt][l - P^], which is the 
probability that N(t) is even and no events occur in [t, t + dt). The lengths 
-dt X < df are uniformly distributed with density X/2. The mean and variance 
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of dT(t) art! given by 


E[dT(t) ] = - dt[l - Xdt] P q + dt[l - Xdt][l - P Q ] 


= dt[P - P ] + 0(dt) 
e o 

and 

var[dT(t) ] « 0(dt) 


( 7 - 7 ) 


where P = 1 - P 
e 


11m 

dt-K) 


0(dt)/dt = 


o 

0. 


is the probability that N(t) is even and 0(dt) satisfies 


Considering the characteristic functional for y 


( to ; d ) = E^e 


where E{*} is the expected value operator, the partial differential equation 
describing the evolution of f(y;t) with time can be derived as in Appendix C. 
The resulting equation 


ft f(Y5t;) = ?y te " 2U " f(Y;t) 


( 7 - 8 ) 


has boundary conditions 


lim 

t-»-0+ 


f ( y; t) - A(o) 


( 7 - 9 ) 


which ensures that all paths start at y Q = 0 and 


f(y=t;t) = f (y = -t;t) + Xe ^ 


( 7 - 10 ) 
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which Insures that the probability mass for values of -t < y < t integrates 

, - -Xt . 

to 1 - e , i.e. , 



f(y;t> dy 



(7-11) 


The rest of the probability mass (e ^ t ) is atomic mass located at y m t, rep- 
resenting the probability that no Poisson event occurs. 

Equations (7-8) through (7-11) summarize the evolution of the probability 
density function for the cross-path deviation y in a Poisson-driven dynamical 
model for aircraft trajectories. 


7.3 Extensions of the Poisson-Driven Dynamical Model 

In order to apply this type of procedure to the analysis of the un- 
controlled terminal area environment it would be necessary to increase 
dimensionality. Cross-path deviations could then be represented in a 
statistical description dependent fundamentally on a Poisson rate parameter. 
Validation can be accomplished through comparison of these descriptions with 
that obtained from the real data. Path structure definition will then be a 
prime factor in the time evolution of these path deviation density functions. 
Once path variations are characterized techniques such as those discussed in 
previous chapters can be used to evaluate hazard. 
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CHAPTER 8 


SUMMARY AND RECOMMENDATIONS 

This report describes results of a preliminary effort in a study 
designed to evaluate prescribed procedures used in the general aviation 
uncontrolled terminal airspace. Fiocedure evaluation is primarily from 
the standpoint of mid-air collin ujn hazard in the "see-and-avoid" environ- 
ment. A generalized expression for mid-air collision (MAC) probability is 
developed. This analysis indicates that traffic pattern design should 
minimize encounter rates and maximize visible avoidable encounters under 
the constraint of pilot visibility inherent in aircraft design. 

Extensive data have been collected by NASA Wallops personnel in several 
uncontrolled environments subject to current prescribed procedure. Some 
results of analysis of these data have indlcate u that pilot adherence to 
present procedure does vary considerably, ludicet ivs are that this varia- 
tion from procedure is based to some extent on individual attempts to maxi- 
mize the ability to "see-and-avoid" within the te/minal airspace. These 
data -re to provide a basis for evaluating procedure in the uncontrolled 
environment . 

Two basic approaches' are described to the evaluation of procedure. 

One involves a deterministic model used to generate pseudo-random paths 
consistent with a given procedure. The ability to define these paths in a 
space-time coordinate system is demonstrated. Furthermore, a method of 
examining the mid-air collision hazard conditioned on a deterministic path 
structure is developed. An approach to the extension of this method to 
random paths is also discussed. A preliminary deterministic model is used 
to demonstrate how validation can be accomplished using real data analysis 
results. Terminal area density of aircraft for the approach phase is 
varied and a hazard measure is calculated to provide information relative 
to the evaluation of a path structure. 

A second procedure developed involves the use of a statistical model. 
This is discussed in terms of a dynamical model which can generate path 
deviations in the form of a stochastic process generating function. 
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Expressions for the probability density function of these deviations are 
derived. This method can be used in conjunction with the developed pro- 
cedures for hazard measures in a manner similar to that used with the 
deterministic model. 

For the continuation of this effort, several recommendations are 
offered based primarily on results of studies thus far. 

(1) NASA Wallops data should be made available in a format and 
organizational structure which would facilitate extensive additional 
analyses. As an example the major analysis done so far includes distri- 
butions of cross-path deviations at several geometric planes within the 
approach pattern. Since no time correlation is available it is of interest 
to do analysis of hazard probability with some assumed time structure. 

This analysis could be the form of statistics such as encounter rates, 
duration of encounters, visibility of hazardous aircraft, etc. 

(2) The available data can provide baseline e^luation of procedures 
to the extent that traffic volumes vary considerably and hazard evaluation 
can be made for a variety of procedures. In the current study time correla- 
tion of the data from individual aircraft has been assumed and is considered 
to be of definite benefit for a complete assessment of the uncontrolled 
airspace. Other information concerning aircraft simultaneously in the 
airspace can definitely complement the analysis, e.g., aircraft mode (landing, 
departure, etc.), and whether path deviations from procedure are pilot 
preference or are in fact the result of avoidance maneuvers. Therefore* 
consideration should be given to acquisition of additional data. 

(3) The deterministic model should be extended to include the total 
terminal airspace so that a complete representation of current procedures 

is available. This procedure should be extended to evaluate various approach 
path structures. Consideration should be given to include in the model the 
decision algorithms necessary to include avoidance maneuvers based on the 
see-and-avoid concept. This can provide more realistic evaluation of 
procedure. Some information is available and more can be obtained con- 
cerning pilot workload, what portion of the time the pilot spends looking 
for other aircraft, and how he spends his look time. Ideas presented in 
this preliminary study indicate that path characteristics can definitely 
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enhance the see-and-avoid concept and should be influenced by this consid- 
eration. For example, most efficient use of the pilot's look time can be 
directly linked to the visibility of hazardous aircraft and the magnitude 
of the search region of potential hazardous aircraft. This is tied to the 
approach path structure and a measure of pilot adherence to the structure. 

It seems reasonable to expect that pilot adherence to procedure will be 
directly proportional to the value of the procedure from the standpoint of 
protection against potentially hazardous situations. 

(4) The statistical modeling should be pursued further to provide a 
viable method for evaluating procedure from a hazard standpoint. Initially, 
this would involve extension of the analytical results to a multidimensional 
situation to more accurately represent the real environment in an uncontrolled 
airspace. This can provide a useful too] to evaluate any airspace environ- 
ment with any conceivable pattern structure. Estimates of pilot adherence to 
prescribed procedure can be input in terms of statistical descriptions. In 
fact limits on adherence to proposed procedure could be a parameter used to 
qualify a hazard measure for any particular postulated and evaluated 
procedure. 
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APPENDIX A 


FLIGHT PATH C NERATION FOR SLOW AIRCRAFT 

This appendix describes the algorithms used to generate aircraft 
tracks in the uncontrolled airport environment. Of twelve possible 30® 
entry sectors defined, only two sectors are included In the initial model 
and described here. 

In this development of a geometric model for generating flight paths 
(tracks) for slow aircraft, several assumptions have been made concerning 
the aircraft and their path characteristics. These assumptions are made pri- 
marily to preserve simplicity in the model but are also based on examination 
of data. The assumptions are: (1) the runway is oriented in a north-south 

direction with the threshold (touch down point) at the south end, (2) the 
threshold is defined as the origin for a cartesian coordinate system with 
angles measured positively from the north-south line (runway centerline) in 
a clockwise direction, (3) aircraft tracks initiate at a point which is 3 
statute miles (4.8 km) from the runway threshold, (4) all winds (if any) are 
from the north (0°), (5) all turns are smooth with a turn rate specified at 
the beginning of a run (up to 30° per position update), (6) all random vari- 
ables are considered to be uniformly distributed, and (7) slow aircraft are 
those which have a terminal airspeed less than about 132 mph (195 fps) and an 
approach airspeed ’.ess than about 70 mph (103 fps) (these two speeds are used 
in the model) . 

A.l Geometric Model for Sector 1 

The Sectoi ] is defined as that 30° arc between 0® and 30° of a 3-mile 
radius circle centered at touchdown. Angular measurement is positive-clock- 
wise from the N-S line. Tracks originate with a uniform distribution between 
0° and 30°. The algorithm functions as follows: 

(1) Upon entering the sector an aircraft is aspigned an entry point 
and an initial heading, 6, where 

f(6) = .1 ; 235° < fl < 245° . 

Note: f(-) indicates probability density function. 
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Figure A-l. Geometry Associated with Sector 1 Initial Entry Segment, 
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(2) The initial leg length, is related to a randomly selected angle, 
a, which is measured counter-clockwise from the N-S line with vertex at touch- 
down. Then 


f (a) - .08333 ; - 22° < a < - 10° 

is used to select a which determines the end of the initial leg, where a 

turn to the downwind leg is initiated. This geometry is illustrated in 
Figure A-l. 

For any selected a the length £,^ iB calculated according to Figure A-2 
as follows. Angle 8 is the actual sector angle at which the entry point is 
located. With initial heading 0, then 

© x = 0 - (180° + 0) 


The angle i|i 


* - (8 - a + 0j) - (0 - a - 180°) 


Therefore the leg length 


l 


1 


n 

COSl|( 


where 


n - r^ sin (0-a) 


yielding 


r sin (B-a) 
*1 - cos (0-a) 


(A-l) 


Thus once a, 8, and 0 are generated the Initial leg length Z ^ may be calculated 
from (A-l) . The distance R given in Figure A-2 which is the distance from the 
terminus of to touchdown is defined as the range of the aircraft. 
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(3.) Upon completion of the Initial leg a new heading HDGDW Is determined 
as the downwind leg heading. A random value y drawn according to 

f(Y) - ; - 10° < Y < 17° 


is used to form 


HDGDW » 180° + y . 

(4.) The downwind heading is flown until y - DWN2 where DWN2 is selected 
from a uniformly distributed range of values - 2500 ft <_ DWN2 ^ - 4000 ft 
which terminates the downwind leg south of the runway, 

(5.) Upon completion of the downwind leg the aircraft Initiates a turn 
until a heading of 90° is reached. At this point the distance to the runway 
centerline is calculated and the difference between this distance and the 
radius of turn is the length of the base leg. The radius of turn is deter- 
mined by the values of the rate of turn and the aircraft velocity. Rate of 
turn is dependent upon the initial parameter setting at the beginning of a 
run and the time between position updates. 

(6.) Upon completion of the base leg a turn la initiated to bring the 
aircraft to a heading of 0° along final approach. 

Figure A-3 illustrates the segments of flight for an entry into Sector 
!. Figures A-4 and A-5 illustrate two typical computer generated Sector 1 
entry tracks. Figure A-6 provides a flowchart of the Sector 1 entry algorithm. 

A. 2 Geometric Model for Sector 11 

Sector 11 is defined as the 30-degree arc between 300° and 330° of a 
three-mile radius circle centered at runway touchdown. Angular measurements 
are clockwise from the N-S line (extended centerline of runway). Upon 
examination of the tracks from Sector 11 in the Wallops data base two basic 
types of tracks were apparent: (1) one type consisted of an initial entry 

leg which terminated close enough to the runway to permit direct entry into 
the downwind leg and (2) the other type consisted of an initial leg which 
ended too far from the runway to start a downwind so that a second leg inward 
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Computer Oenerated Track for Sector 1 

i i-.uro A' 4 . Kxamplc’ o" Computer r.cnornLod Arrival Aircraft Track with 
l'.ntry into Sector ] . 
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Computer Generated Track for Sector 1 

ip, ure A- 5. l'xninplc*. of Computer Generated Arrival Aircraft Truck v;ith 
Kntry into Sector 1. 
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Fig. A-6. Program to generate track of slow aircraft arriving from 
Sector 1. 




was necessary to get to the downwind leg. Examples from the actual data 
are shown in Figure A-7. A definition of these two types is shown in 
Figure A-8. 

The simulation model provides for generation of these types of slow 
aircraft tracks entering Sector 11. The description of the algorithm is 
as follows: 

(1.) Coordinates of entry (XAPT, YAPT) are selected at random assuming 
an equally likely distribution been 300° and 330° of the angle AANGL defined 
in Figure A-9. Thus 


XAPT = - RADIUS * SIN (360 - AANGL) 
YAPT = - RADIUS * COS (360 - AANGL) . 


(2.) The initial heading (HDGIN) is uniformly distributed between 118° 
and 180°. It is necessary to have the initial leg length depend on AANGL 
so that the length corresponding to AANGL - 330° will not cause the path to be too 
close to the runway centerline or overshoot it. The entering aircraft must 
always be able to enter a downwind leg. At AANGL = 330° the initial leg 
length is randomly chosen from an even distribution of lengths between 2500 
and 5000 feet and at AANGL = 300° this range is between 5000 and 10000 feet. 

Each leg length limit varies linearly with AANGL as 


LL (lower limit) 


2500 + (330 - AANGL) * 


5000-2500 

330-300 


and 


UL (upper limit) = 5000 + (330 - AANGL) 


* 10000-5000 
330-300 


and is illustrated In Figure A-10. 

(3.) The coordinates at the end of the initial entry leg are calculated as 

SBDWT = XAPT + DISIN * SIN(180-HDGIN) 

YBDWT “ YAPT - DISIN * COS(180-HDGIN) 


as per Figure A-ll. 
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Figure A-7. Actual Tracks of Aircraft Arriving from Sector 11. 
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Runway Centerline 


Touchdown 


(XBFNT, YBFNT) 


(b) TYPE II 


Figure A-8. 


Two Types of Sector 1 1 Entry Approach Paths Modelled. 




AANGL 


Figure A-9. Sector H Entry Point Geometry. 


Leg Length 



Figure A-10. Sector 11 Initial Leg Length Distribution vs AANGL (defined 
in Fig. A-9) . 
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Figure A-li. Sector 1.1 Initial Entry Leg Terminus (XBDWT, YBDWT) . 
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Figure A-12. Transition Regions for Sector U Entry. 



(4.) If the aircraft is within 7500 ft of runway centerline, It la in a 
suitable position to proceed downwind. If not, it la made to proceed further 
inward on a heading of 135° (fixed) until it is within 5000 ft; then, it 
proceeds downwind. This latter situation occurs if lower AANGL values have 
higher HDGIN values (see Figure A-12) . 

(5.) Turns (if any) after initial leg are made at terminal velocity 
(VELHI). The TURNS subroutine accepts old and new headings, velocity, and 
coordinates at beginning of turn; it returns coordinates after turn. 

When the initial leg terminates at a distance from the runway centerline 
greater than 7500 ft, a turn is made to get to a new heading of HDG2L * 135°. 
Downwind leg entry is delayed with new coordinates-before-downwlnd-turn 
calculated as XBDWT * - 5000 ft (see subparagraph 4 above) and YBDWT ■ YAFTT 
+ (5000 + XAFFT) COTAN (180 - HDG2L) . The program then continues with the 
following logic. 

Once the termination of the initial entry leg is within 7500 ft of the 
runway centerline, a downwind heading is generated by sampling between 175° 
and 185° (see Figure A-12) . 

(6.) The downwind leg terminates between 2500 and 5000 ft below the 
threshold. XBBST, YBBST are coordinates of terminus of downwind leg where 
turn onto base leg begins (see Figure A-12) . 

(7.) A base heading is chosen between 085 and 095 from a uniform 
distribution (HBA3E) (see Figure A-12) . 

(8.) The base leg length is determined from the point of termination 
of the downwind leg, the base leg heading (HBASE) and the distance required 
to make the turn to final. Figure A-13 illustrates the geometry involved in 
determining the distance YBENT and YBENT yielding the base leg distance 

DIST = { [ XBENT - X(NN)] 2 + [YBENT - Y(NN)] 2 } 1 ^ 2 

where X(NN), Y(NN) is the beginning of the base leg. From Figure A-13 
XBENT and YBENT can be determined as follows. 

X(NN) - XBBST + R + R sin(90-HBASE) 

Y(NN) = YBBST - [R - R cos(90-HBASE) 
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Figure A- 13. Geometry for BASE LEG and FINAL APPROACH TURN for Sector !3 
Entry . 

XBENT - - [R - R sin(90-HBASE)] 

YBENT « Y(NN) + [XBENT - X(NN)] tan(90-RBASE) (A-l) 

where R Is the turn radius and (9C-HBASE) * 8 as shown in Figure A-13. In 
the model these calculations are obtained using the subroutine TURNS. The 
downwind leg heading, the desired base leg heading, velocity and position 
(XBBST, YBBST) are input to yield [X(NN) ,Y(NN) ] . Then using TURNS with the 
same input, except replacing desired base leg heading with the final approach 
heading, an output position [X(N),Y(N)] is used to find XBENT as 

XBENT - X(NN) - X(N) . 

Once XBENT is found, YBENT can be found using eq. (A-l). 

(9.) The turn to the final approach heading (0°) is made at the same 
slower velocity at which the base leg is flown. Subroutine TURNS is agait 
used. 

Figures A-14, A-15 and A-16 provide plots of typical approach tracks 
generated by the computer algorithm. 

Figure A-17 gives a flowchart of the algorithm. 
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Figure A- Jo. F.xnnip I v of Guntputor (’eneraletl Arrival Aircraft Irack with 
r.nlry into Hoc Lor 11. 


105 



I'ij'.iitv A-l r j. r.v;imp1o of Computer Generated Aircraft Track with 
i -it t r y [ nt o Sector 11 • 


106 



Figure A-16. Example of Computer denerated Aircraft Track with Entry 
Into [Sector 11, 
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Hu. A-17, Program to generate track of slow aircraft arriving from 
Sector 11. 
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APPENDIX B 


AIR TRAFFIC SIMULATOR 

The air traffic simulator program to study the uncontrolled airport 
environment currently consists of five major programs as outlined in 
Figure B-l. A brief description of each program along with flowcharts 
in Figure B-3 follows. 
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Fig, B-I. General flowchart of air traffic programs, 
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Fiq. B-l. Continued. 
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Fig. B-l. Concluded. 




PROB Is a program which incorporates two logical units; the Sector 
Selector and the Output Unit. The Sector Selector unit of PROB is executed 
if the number of paths generated has not exceeded the desired total. The 
function of this unit is to assign a value to the time parameter, TAU, and 
to select a sector based on a predetermined probability of selection. 

TAU, the inter-arrival time between aircraft 1-1 and aircraft 1, is 
computed by the following equation 

TAU - - j £n(x) , 

X « number of arrivals/time period 
x = uniform random number. 

The Output unit of PROB is executed after all paths are generated. It 
calculates the final mean and standard deviation for each desired distribution 
and prints the results along with the histograms. 

DISTR incorporates the Statistics Unit. Execution of this program 
occurs after each flight path has been generated. Its function is to 
determine by interpolation where and when the aircraft crossed a given 
geometric plane. Figure B-2 illustrates this for a given set of points. 


where 

and 


• Pol.ir (0 ,y , t) where y and t are found 

Point 1 q t>y erpolating between point .1 and 

(xi>y^>ti) " s pole ’. 

1 — ± — i x = 0 

•n 

Point 2 

(x 2 , y 2» t; 2 ) 


Figure B-2. Sequential Positions of Aircraft with Respect to Plane 
of Interest at x = 0. 
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If the flight path is the very first one to ba generated than it la 
not Included in the statistics but serves as a basiB from which to obtain 
data for the next flight. Parameters TARG and IDSK are used to determine 
which statistics are to ba maintained. 

WSL01 is a Path Generator program (see Appendix A). It handles tha 
aircraft entering Sector 1 (0 to 30 degrees of N-S line). Paths are 
generated in either straight line or curved segments. WSLOl determines 
the end points of each segment of the flight path. Intermediate points 
along the segment are obtained by calling subroutine FLY2. 

WSLll is a Path Generating program. It handles air 'raft entering 
Sector 11 (300-330 degrees). Paths are generated in segments as in WSL01, 
and the subroutine FLY2 calculates intermediate points along the segment. 

A detailed description of the trajectory algorithms appears in Appendix A. 

FLY2: This subroutine called from any Path Generating program calculates 

intermediate points along a segment. Provision is made for partial time 
intervals and for straight or curved segments. 

The computer model has the following limitations: 

1) It currently handles aircraft in only two of twelve sectors. 

2) Aircraft altitude is not considered. 

3) There are no provisions for aircraft flying by or taking off. 
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Fig. B-3. Individual flowcharts for air traffic programs. 
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Fig. B-3. Continued. 
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Fig. B-3. 


Continued. 
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Fig. B-3. Continued. 
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Fig. ii-i. Continued. 
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Fig, H-3. Continued. 
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Fig, 13-3. Concluded. 
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APPENDIX C 


DIFFERENTIAL EQUATION FOR CROSS-PATH DEVIATIONS IN POISSON-DRIVEN 

DYNAMICAL SYSTEM 

Consider the random differential equation describing cross-path 
deviation y at time t 

d + S dt - dT(t) (C-l) 

y t y t 

where 

y t . e y t + / . - S(t l- a) (-1) N<0) do 

1 4 

where (3 is a damping term. Here y is an integrated version of the random 
telegraph wave (ref. 11) described by 

T(t) = (-l) N(t) 

where N(t) is a Poisson counting process defined as a control (+1) between 
event times of the Poisson process. 

Let f(y;t) denote the first order density function for y^. The incre- 
ment dT(t) can be defined as 

t+dt dt 

dT(t) - J' (-1) N(0) do = j (-l) N(t+U ^ du 
t o 

- - (dt - 2a) (-1)^'*"^ , 0 < a < dt 

where a is an event time. Consider the characteristic functional for y^ 

1 1 
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where E{ • } is the expected value operator. This can be expressed as 


>;t) - f e j(11Y f (y;t) 


The partial derivative with respect to t is 


J dy It f (Yit) 


This result can also be obtained by examining 


— - “" 0 5(»;t + At) - J(m| 


lin -i T ia7 t + 4y t J“y t 
' ;;!n «t> 1 E|e £ ‘ - e C J 


lim „,.-l J J“ y t / 3 “ 4y t . ^1 

At-*) <4t> E L _1 /J 


The expectation can be written as a conditional expectation 


a$(ai;t) 

3t 


- j „ J “ y t lin, ,.,.,-1 „ J “ 4y t ,, 

*y. ) e ^ (aa) E -i) y. 



where E 
expectaSl 


denotes expectation under the distribution of y 
on in eq. (C-3) can be written 


f 


The latter 


E 






+ 



From the equation for dy £ given in eq. (C-l) , 


' J 

E (e 


juAy. 


-1) y 


t. = 


By*. + p c 


- V 


At + 0 (At) 


Therefore 


- If - 1 - By t | e j Vt jm[P e - P o - eyjj • 

A final simplification follows by observing that the expected value of the 
random telegraph wave satisfies 

E[T(t)] - P - P 
e o 

It is easily shown that (ref. 12) 


E[T(t) ] = e “ 2At . 


Thus we can write 


■/ 




t , r -2Xt 

jui [e - (Jy 


.} 


dye ju [e 2Xt - 0 y 3 f(y;t) 


CC-4) 
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Equating eqs. (C— 2) and (C-4) , yields 



It 





or, integrating the right-hand-side by parts 


■fjr f (y;t) ** —• [e“ 2Xt - 3y] f(y;t) 


(C-5) 


This partial differential equation describes the evolution of f(y;t) with 
time, 
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